Radiative capture and electromagnetic dissociation involving loosely bound nuclei: 

the 8 B example 
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Electromagnetic processes in loosely bound nuclei are investigated using an analytical model. In 
particular, electromagnetic dissociation of 8 B is studied and the results of our analytical model are 
compared to numerical calculations based on a three-body picture of the 8 B bound state. The 
calculation of energy spectra is shown to be strongly model dependent. This is demonstrated by 
investigating the sensitivity to the rms intercluster distance, the few-body behavior, and the effects 
of final state interaction. In contrast, the fraction of the energy spectrum which can be attributed 
to El transitions is found to be almost model independent at small relative energies. This finding 
is of great importance for astrophysical applications as it provides us with a new tool to extract 
the El component from measured energy spectra. An additional, and independent, method is also 
proposed as it is demonstrated how two sets of experimental data, obtained with different beam 
energy and/or minimum impact parameter, can be used to extract the El component. 

PACS numbers: 21.60.Gx, 25.60.Dz, 25.70.De, 27.20.+n 



I. INTRODUCTION 

The properties of loosely bound nuclei have been stud- 
ied in nuclear physics for a number of years. In par- 
ticular some electromagnetic processes, such as certain 
charged particle capture reactions, are very interesting 
in themselves as they are of vital importance in astro- 
physical scenarios. Unfortunately, at stellar energies, the 
cross sections for these reactions are very small due to the 
Coulomb barrier and direct measurements are therefore 
extremely difficult. Instead one has to rely on theoret- 
ical extrapolations from experimentally accessible ener- 
gies down to stellar ones. An alternative, indirect method 
to investigate radiative capture reactions, is to study 
electromagnetic dissociation (EMD) on heavy targets [1]. 
This technique gives an enormous increase in yield due to 
the huge amount of virtual photons produced by the high- 
Z target, the more favorable phase space factor, and the 
possibility to use thicker targets. In principle it should 
therefore be possible to measure the cross section at very 
low relative energies. There are, however, also disad- 
vantages with the indirect method; the most important 
being the admixture of 7-transitions with different mul- 
tipolarities whereas the direct capture process, in most 
cases, should be completely dominated by El transitions. 
Since the cross sections for radiative capture and photo- 
dissociation are related via detailed balance for each sep- 
arate multipolc, it is necessary to have knowledge of the 
strengths of different multipole transitions in the EMD 
reaction. 

The problem of extracting the El contribution from 
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a measured EMD energy spectrum remains a challenge 
to the nuclear physics community. In Ref. [2] it was pro- 
posed to study the angular or momentum distributions of 
the breakup fragments. The idea was to employ the fact 
that interference terms, between El and E2 excitation 
amplitudes, will produce an asymmetry in these distri- 
butions. This method was, e.g., used in the analysis of 8 B 
EMD [3] where the E2 excitation amplitude, calculated 
within first-order perturbation theory, was renormalized 
in order to reproduce the asymmetry of the measured 
7 Be longitudinal momentum distribution. However, as 
was already noted in Ref. [2], the asymmetry due to El- 
E2 interference strongly depends on the final state inter- 
action (FSI) between the breakup fragments; or in other 
words, on the structure of the continuum up to relatively 
large energies. Moreover, terms which contribute to the 
asymmetry do not contribute to the integrated cross sec- 
tion from which the astrophysical iS-factor is extracted. 
Finally, for low beam energies, higher-order dynamical ef- 
fects will lead to a reduction of the asymmetry [4] . There- 
fore we conclude that, if we are interested in astrophysical 
applications of EMD experiments, it is desirable to look 
for more stable, and less model dependent, characteristics 
than the asymmetries of angular and momentum distri- 
butions. In this paper we will present two novel methods 
to extract the E1/E2 ratio from EMD experiments. 

We will use an analytical approach based on a two- 
cluster picture of the nucleus, but the effects of many- 
body structure will also be included. Our approach 
will be general in the sense that both neutron-rich and 
proton-rich systems can be studied. This model was first 
presented in [5] while similar approaches also exist for 
one-neutron [6, 7] and two-neutron [8-10] halo nuclei. 
Although advanced numerical investigations are readily 
performed utilizing present day computer power, an ana- 
lytical approach might have an advantage when exploring 
general physics features, and the sensitivity to different 
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model assumptions. 

The physics case which will be investigated throughout 
this paper is the 8 B nucleus. The interest in 8 B stems 
from its key role in the production of high-energy solar 
neutrinos. It is well-known that the probability for the re- 
action 7 Be(p, 7) 8 B at solar energies strongly depends on 
the structure of 8 B and, in particular, on the asymptotics 
of the valence proton wave function (WF). This reaction 
has been studied indirectly through EMD, using a ra- 
dioactive 8 B beam impinging on a heavy target [3, 11-13]. 
Note that the question of E2 contributions to the exper- 
imental spectra was addressed differently in all these in- 
vestigations. We should also mention the recent progress 
in radiative capture measurements [14], where the cross 
section has been measured at energies around 200 keV 
with an accuracy of « 15 percent. Nevertheless, in all 
cases theoretical models are needed to extrapolate the 
measured cross sections down to solar energies. Theoret- 
ical studies of the low-energy behavior of the astrophys- 
ical S'17-factor has been presented by many authors, see 
e.g. [15-21]. 

The structure of this paper is the following: Section II 
contains a summary of the theoretical framework that 
will be used in the calculations. In Sec. Ill our analytical 
model WFs are presented and discussed in quite some de- 
tail. Finally, in Sec. IV we discuss the model dependence 
of calculated EMD energy spectra and propose two new 
methods to extract information from EMD experiments. 



II. THEORETICAL FRAMEWORK 

Our starting point for calculating electromagnetic 
cross sections will be the EA strength function for a tran- 
sition from a bound state (total spin Jo) to a continuum 
state with energy E 



dg(EA) 
dE 



2J + 1 



*E/ dr 3 \(j\\M(E\)\\0)\ 2 S(E 3 -E), (1) 
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where drj is the phase space element for final states, 
M(EX, fi) is the electric multipole operator and |0), \j) 
are the bound and continuum states in the center of mass 
system. 

We will consider loosely bound systems of two clus- 
ters (c + x) and, in particular, we will study transi- 
tions to the low-energy continuum in which excitations 
are manifested as relative motion between the clusters 
E = h 2 k 2 /2fi cx , where [i cx is the reduced mass of the 
two-body system. Introducing the intercluster distance 
r, the corresponding cluster EA operator (operating only 
on the relative motion of clusters) is 



M(E\,») = eZ(\)r x Y Xft (r), 



(2) 



where we have also introduced the effective multipole 
charge Z{\) = (i x x (Z x /m x + (-l) x Z c /m x ). 



The strength function is the key to study several 
reactions. For example, the cross section for photo- 
dissociation A(j, x)c is given by 



EA 



(E) = 



(2tt) 3 (A + 1) (E^- 1 dB(EA) 
A[(2A + 1)!!] 2 



dE 



(3) 



where the photon energy E y = E + E is larger than 
the binding energy Eq. From this formula the inverse 
radiative capture reaction c(x, j)A can be studied using 
detailed balance 



(E) 



\hck 



2(2 J A + 1) 



(2J C + 1)(2J X + 1) 



< A (£), (4) 



where Jj is the spin of particle i. Note that the probabil- 
ity for direct capture of charged particles is dramatically 
reduced at low energies due to the Coulomb barrier in 
the c + x channel. The cross section is therefore usually 
factorized into the Gamow penetration factor and the 
5-factor 



a IC {E) = 



-2-rrn(k) 



E 



-S(E), 



(5) 



where r](k) = Z c Z x e 2 [i cx /h 2 k is the Sommerfeld param- 
eter. The dominant part of the energy behavior is car- 
ried by the Gamow penetration factor while, e.g., nuclear 
structure information is incorporated into the S'-factor. 

Finally, we will consider the process of EMD on a 
high-Z target. Using first-order perturbation theory, and 
the method of virtual quanta [22, 23], the energy spec- 
trum can be written as a sum over multipole (nX) photo- 
dissociation cross sections multiplied by the correspond- 
ing spectra of virtual photons n^^E^) 

d(7EMD n E\(E 7 ) eA/t-hn 

= 2^— — a -i ( E > 



dE 



E~ 

EA 7 



+ 



~E, 

MA 



<\E). (6) 



Note that, except in the vicinity of corresponding 
resonances, MA transitions arc usually strongly sup- 
pressed [23]. Therefore, we will not study them in this 
work. 



III. ANALYTICAL MODEL 

A. Model wave functions 

A straightforward calculation of the electric multipole 
matrix element for a direct transition between a loosely 
bound state and a non-resonant continuum state, shows 
that the radial integrand rises to a maximum value at a 
radius which is, in most cases, many times the nuclear ra- 
dius. Thus, these processes will mainly probe the surface 
structure of the nucleus. Furthermore, "loosely bound" 



3 



implies that the nucleus will exhibit a large degree of 
clusterization and that the relative motion WF between 
the core and the valence nucleon will have an extended 
tail. 

The final, continuum state will contain both Coulomb 
and nuclear distortions. For low continuum energies, and 
when the binding energy of the initial state is small, the 
nuclear distortions can be neglected in a first approxima- 
tion. Therefore, we will only consider a pure Coulomb 
continuum in our analytical model, i.e., all nuclear phase 
shifts will be put equal to zero. The effects of nuclear 
distortions in 8 B EMD will be considered in Sec. IV A. 
Thus, a continuum state, with relative orbital momentum 
I between the clusters, will be described by a normalized, 
regular Coulomb function 



<j>l{k,r) 



7T k 



Fi(k,r), 



(7) 



where 



F l (k,r) = C l (ky kr (kr) l+1 

x 1 F 1 (l + l + it](k);2l + 2;-2ikr), (8) 

and g i is the Coulomb phase, rj(k) is the Sommerfeld 
parameter, iFi(a;b;z) is the confluent hypergeometric 
function [24], and 

d(k) = 2 l e-^' 2 \T(l + 1 + i V (k))\/(2l + 1)!. (9) 

The reduced matrix element introduced in the defini- 
tion of the strength function, Eq. (1), contains a radial in- 
tegral. With our approximation for the continuum state 
this integral takes the form 



h{k) = / dre- lkr r l+1 
Jo 

x iiq(Z + 1 - ir)(k);2l + 2;2ikr)r x 4>(0,r). (10) 

Here, 0(0, r) is the two-body, relative motion WF de- 
scribing the initial, bound state. At large r, with relative 
orbital momentum / between the clusters, this radial 
function should be proportional to the Whittakcr func- 
tion W_ 71a j a+1 / 2 {2n Q r), see e.g. Ref. [24], where 770 = 
Z c Z x e 2 /j cx /h 2 Ko and E — ti 2 K 2 l /2fi cx is the binding en- 
ergy 

In most studies on loosely bound systems, the Whit- 
taker function has been used to describe the bound state 
for all r. However, the Whittaker function behaves as 
r~ l ° in the limit r — -> 0, and therefore this approxima- 
tion is only motivated if the transition matrix element 
is dominated by contributions from very large r. This 
is the case for reactions at very small energies; while for 
real experimental energies (E > 100 keV), the WF of 
the bound state should be constructed in a more realistic 
way. Our idea is therefore to introduce a model function 
that describes the bound state (c + x) WF accurately for 
all distances. This can be achieved by considering the 



behavior at small and large r. We have already pointed 
out that the WF should be described by a Whittaker 
function at large r. Furthermore, the expected r — > 
behavior for a two-body system consisting of point-like 
particles is r lo+1 . Both asymptotics are fulfilled using the 
following model function 



<^ cxact "(0,r) 



1 



--W- Vo j 0+1/2 (2K Q r) 



x (1 - e-^ r f 0+ \ (11) 



where is the normalization constant and 7 denotes 
the parameters {ko, k\, 770}. The parameters kq and 770 
are defined by the binding energy, charges, and masses, 
while Ki can be fitted to give the correct distance between 
particles c and x (or the correct size of the system). Us- 
ing this WF, and solving the integral (10) numerically, 
it is possible to get very good estimates for the electro- 
magnetic reaction cross sections. We should also mention 
that in the limit fc^Owe are actually able to solve the 
radial integral (10) analytically, except for the normal- 
ization constant N^. 

However, we are searching for a completely analytical 
model which will also enable us to incorporate many- 
body effects. Our model function has to be modified 
accordingly. First, we note that the asymptotic form of 
the Whittakcr function as r — > 00 is 



W_ VOtlo+1/2 (2K r) 



{2n r)^ 

(Vo - jo) + fog - ll) 
2n r 



• (12) 



Secondly, for two-body systems in which the clusters have 
an internal structure, the centrifugal barrier is effectively 
larger and the WF should behave as r n (where n > Iq + 1) 
as r — > 0. 

Motivated by this, we put forward the following model 
function 



y(0,r) = 



1 



(13) 



with norm 



2p 



N ~< = E m(-i) m (2«;o+m Kl ) 2 ^- 1 r(i-2^), (14) 



m=0 

where 7 denotes the parameters {k , m, ?7o,p}, and p is 
an integer fulfilling p > rj' + 1. The parameter k is 
defined by the binding energy and effective mass. By 
putting rj' = ijo we would ensure to reproduce the tail of 
the WF at very large r. However, the difference between 
an exact Whittaker function and its asymptotic behavior 
(first term of Eq. (12)) remains important for r < 100 fm. 
Therefore, r]' and K\ are used as free parameters in a fit 
to the "exact" WF (11) in the interval of interest. In this 
way i]' will be an effective Sommerfeld parameter while 
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K\ will still mainly be connected with the size. Note 
that if r;o > lo, then the second term in Eq. (12) will 
be negative and consequently we will find that t]q < r) . 
Finally, the integer p is fixed by the small r behavior. For 
a pure two-body system we will use p =< r]' + l + 1 > 
(where < x > is the closest integer to x), while we can 
take many-body effects into account by putting p =< 
r)' + n>. 

With this model WF it is possible to solve the inte- 
gral (10) exactly 



Ii tl {k) 



1 



E 

"7 m =0 



m 



i-iy 



x (mm +n + ik)-V +2+x -^T{l + 2 + A - r]' ) 
x 2 Fi + 2 + A - r]' ; I + 1 - ir](k); 



2/ + 2; 



2ik 



mKi + kq + ik 



(15) 



Many-body nuclear structure can further be taken into 
account by considering the possibility that the bound 
state WF contains several different two-body components 



0(0, r) = ^a^ 7i (0,r), 



(16) 



which can be seen as two-body projections of the many- 
body WF. Note that pure many-body components will 
not contribute to two-body breakup and, as a result, we 
will have JA of < 1. Note also that the threshold for 
two-body breakup will be higher for components where 
one (or both) of the clusters is excited. Therefore, we de- 
fine the continuum strength function separately for each 
component. Finally, we arrive at an analytical formula 
for the strength function of component i 



dB(EA) 



e 2 Z 2 (\)[i cx 2A + 1 



dE . H 2 2ir 2 

x^a 2 k 2l+1 tf{k){l Q\Q\lQ) 2 \Ii^(k)\ 2 . (17) 
i 

B. One-neutron halo systems 

The special case where Z x = 0, i.e., a one- neutron halo 
system will lead to several simplifications. First of all we 
will have 770 = and the Whittakcr function in Eq. (11) 
will transform into a modified, spherical Bessel function 



Wo,z 0+ i/2(2K r) 



K la+1/2 {n r). (18) 



Furthermore, the continuum solution will reduce from a 
Coulomb function (7) to the corresponding component of 
a plane wave 



(19) 



0, (ln) (fc,r) = \ -i l rji(kr), 



where ji(x) is a spherical Bessel function. In this case, 
the integral (10) can be solved exactly. Consider, e.g., a 
node-less s state, for which our model WF would read 



,(ln,ls) tr . \ 

<i>\ (0,r) 



1-e" 



with norm 



(In, Is) 



7 



2k («i + 2/c )(«i + «o) ' 



(20) 



(21) 



The radial integral (10) for this special case will be given 

by 



l^\k) 



1 (20!!(/ + l + A)! 



1 

x E( _1 ) ro ( m/tl + «o)" (i+2+A) 



^ ,' 1 + 2 + X 1 + 3 + X 3 -k 2 
< ~ ! " 1 2 5 2 ; 2 + ; (m«i + « ) 2 



C. Application to 8 B 



(22) 



We will now apply our model to the 8 B nucleus. The 
low-lying 8 B continuum can, with relatively good preci- 
sion, be approximated as a pure Coulomb one. At least 
there are no negative parity states at low excitation ener- 
gies [25-27] and the electromagnetic processes are, in all 
cases we are considering, dominated by El transitions. 
However, we will also show that the influence of a broad 
negative parity state at high excitation energy [26, 27] 
is not negligible at intermediate (> 0.4 MeV) continuum 
energies. 

In a first approximation the bound state of 8 B can be 
treated as a pure two-body ( 7 Be+p) system with binding 
energy Eq = 137 keV and relative orbital momentum 
lo = 1. The single free parameter ki in our "exact" 
model function (11) was then fitted to an rms intercluster 
distance of r rms = 4.57 fm (extracted from Ref. [28]). In 
order to get analytical results, we then introduced the 
model function (13). We put p = 3 and fixed k from the 
binding energy, while the remaining parameters «i and 
r}' were fitted to the behavior of the "exact" model WF, 
see Table I. The resulting El strength function is shown 
as a dotted line in Fig. 1. This analytical approximation 
agrees very well with numerical results obtained keeping 
the "exact" model WF. The error is less than 2% in the 
region of interest. 

However, concerning the structure of the 8 B ground 
state one should keep in mind that the 7 Be core is in 
itself a weakly bound system with an excited l/2~ state 
at 429 keV. The common treatment of 8 B as a pure two- 
body system is therefore questionable. We want to inves- 
tigate what effect the many-body structure of 8 B might 
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FIG. 1: The El strength function of 8 B. Although the to- 
tal strengths are the same in these examples, the shapes are 
very different when treating the 8 B nucleus as a two-body 
(dotted line) or a three-body (solid line) system. The differ- 
ence is mainly due to the stronger centrifugal barrier in the 
three-body case (see, e.g., Ref. [29]) which forces the WF to 
be narrower in coordinate space and thus wider in momen- 
tum/energy space. The parameters of the model WFs can be 
found in Table I. 



TABLE I: Parameters of our model WFs used to describe 
the 8 B ground state. Both models give the same intercluster 
distance, r rms = 4.57 fm. With Eo = 137 keV we get Ko = 
0.076 fm -1 and 770 = 1.595. The excited core component (last 
row) has Eo = 566 keV giving Ko = 0.154 fm -1 and ijo = 
0.786. The relative orbital momentum for all components are 
lo = 1 while / is the channel spin and a 2 is the spectroscopic 
factor. Note that there is no dependence on the particular 
value of the channel spin in the two-body case. 



Model WF 


configuration 


1 


2 

a 


P 


Ki (fm" 1 ) 


Vo/Vo 


two-body 


['Be(3/2-)®p] 


2 


1.00 


3 


0.601 


0.79 


three- body 


[ Y Be(3/2~)®p] 


2 


0.65 


5 


0.702 


0.87 




[ 7 Be(3/2-)®p] 
[ 7 Be(l/2-)®p] 


1 
1 


0.13 
0.16 


5 
5 


0.765 
0.753 


0.86 
1.43 



D. Studies of the corresponding potential 

Using the two-body WF (11), which describes correctly 
the binding energy and the geometry, we are able to re- 
store the corresponding two-body potential. Besides cen- 
trifugal barrier and Coulomb interaction, the potential 
also contains an attractive part which we find can be ap- 
proximated with a high accuracy by one or two Yukawa- 
type potentials. Note that such potentials are widely 
used in few-body nuclear physics. In Fig. 2 the nuclear 
part of the potential which corresponds to our 8 B "exact" 
two-body WF is plotted. As is shown in the figure, the 
potential can be very well described by a Yukawa-type 
potential. In this connection we would also like to point 
out that for the special case where lo = and Z x = 0, 
the WF can be described by (20) and the corresponding 
potential reduces to a Hulthen potential which has exact 
solutions. 



> 




r (fm) 

FIG. 2: The nuclear potential (dotted line) corresponding to 
our 8 B "exact" WF. This potential is very well described by 



a Yukawa potential V n - 
line) . 



z (r) = -114e" 



r /r (MeV) (dashed 



IV. ELECTROMAGNETIC DISSOCIATION 



have on the strength function. For this purpose we utilize 
a recent three-body (a + 3 Hc + p) calculation [28] where 
it was shown that, after projection onto the two-body 
channels, there are three main components (adding up 
to 94% of the total WF, see Table I) and that the rest 
are pure three-body channels. For each of the numeri- 
cal two-body overlap functions we fit our parameters m 
and 7/q. The binding energy, Eo = 137 keV, determines 
ko for the two first components and Eo = 566 keV for 
the third, 7 Be excited state, component. The best fit 
to the small r behavior is obtained with p = 5 which 
reflects the effectively larger centrifugal barrier in the 
three-body case. This centrifugal barrier will push the 
WF away from r = and will, thus, force it to become 
more narrow than the corresponding two-body WF. 



The EMD of loosely bound nuclei, impinging on a high- 
Z target, has been used in nuclear physics for many years 
both in order to investigate nuclear structure and as an 
indirect method to extract information on radiative cap- 
ture reactions. Unfortunately, as for all reaction exper- 
iments, a lot of information is contained in the experi- 
mental results and it is a hard task to disentangle the 
desired part. The transition matrix elements represents 
the probability for an initial state wave function to end 
up in a specific final state after being filtered through 
the reaction mechanism. Naturally, both the structure of 
the initial as well as the final state are important for this 
quantity. In addition, further complications arise if nu- 
clear induced breakup contributes to the measured cross 
section, and/or if the interaction time is long enough for 
higher-order transitions to become important. There- 
fore, in order to minimize interference from higher-order 
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dynamical effects and from nuclear interactions, we will 
be interested in high-energy experiments in which events 
characterized by large impact parameters have been se- 
lected. 

In this section we will demonstrate that the analysis of 
EMD energy spectra from loosely bound nuclei is highly 
model dependent. We will then discuss the important 
issue of how to separate the contributions from specific 
multipoles; a problem which is of great significance when 
extracting information on the inverse, radiative capture 
reaction. 



A. Model dependence of energy spectrum analysis 

Applying our analytical model to study EMD using 
first-order perturbation theory enables us to investigate 
the sensitivity of the energy spectrum to some properties 
of the initial bound state. At small relative energies the 
electromagnetic processes are highly peripheral, which 
means that the interaction mainly probes the external 
part of the bound state WF. As a consequence, the am- 
plitude of the EMD cross section should depend crucially 
on the size of the nuclear system. A larger size implies a 
lower Coulomb barrier, which results in a larger tunnel- 
ing probability, and consequently a larger cross section. 
However, the radii of nuclei far from stability are usually 
extracted from interaction cross section measurements, 
and this procedure has unfortunately proven to be highly 
model dependent. A Glauber- type analysis, assuming a 
uniform density distribution, results in a smaller radius 
as compared to an analysis in which the few-body struc- 
ture is taken into account explicitly [30]. Furthermore, 
the relevant parameter for two-body breakup is in real- 
ity the intercluster distance rather than the total matter 
radius, and the relation between these two quantities is 
also model dependent. In a pure two-body model one of- 
ten assumes that the size of the core is equal to the size 
of the corresponding free nucleus. In contrast, taking 
many-body structure into account will result in polariza- 
tion effects. For example, it was found in Ref. [31] that 
the average distance between 3 He — a is approximately 
10% smaller inside 8 B (studied in a 3 He + a+p picture) 
than in a free Be nucleus ( 3 He + a picture). 

We have investigated the sensitivity of the EMD cross 
section to the rms intercluster distance by performing 
model calculations with a 8 B-like system. A pure two- 
body system with Iq = 1 relative motion (having unity 
spectroscopic factor) was assumed and a binding energy 
of 0.137 MeV was used. The effective Sommerfeld param- 
eter was fitted to the asymptotic behavior of (11), i.e., 
to the Whittaker tail. The size of our model WF (13) 
was then a function of the remaining free parameter n\. 
By varying this parameter we could investigate the sen- 
sitivity to the size and from Fig. 3 it is clear how strong 
it is. Just for illustration: a 15% uncertainty in the rms 
intercluster distance (which, in our model, would corre- 
spond to 5% uncertainty in the 8 B matter radius) results 



in ~ 50% ambiguity of the calculated total cross section. 




FIG. 3: The total EMD cross section as a function of inter- 
cluster distance. The parameter m of our two-body model 
WF was varied in order to change the rms distance between 
the 7 Be core and the valence proton. A beam energy of 
82.7 MeV/nucleon and a minimum impact parameter of 30 fm 
was used. 

Let us now consider the difference between a two- 
body and a three-body approach. As was mentioned 
in Sec. IIIC, the effectively larger centrifugal barrier in 
a three-body system will push the relative motion WF 
away from r = and consequently force it to become 
more narrow than the corresponding two-body WF for 
a given radius. We therefore expect the distribution in 
momentum/energy space to be broader. This effect is 
clearly seen in Fig. 1 where the 8 B El strength functions, 
obtained using our three-body [32] and two-body analyt- 
ical model WFs, are compared. This difference, seen in 
the strength function, should be even more pronounced 
in the energy spectrum since it will be magnified by the 
spectrum of virtual photons. 

In Fig. 4 we compare different calculations of 8 B EMD 
on Pb, including both El and E2 transitions, to the ex- 
perimental data from Davids et al. [3]. This experiment 
is very appealing since the selection of scattering angles 
(8s B < 1.77°, which corresponds to a minimum impact 
parameter of & m ; n = 30 fm) minimizes the contribution 
from nuclear scattering, and the relatively high beam en- 
ergy (82.7 MeV/nucleon) justifies the use of first-order 
perturbation theory. Let us first compare our analyt- 
ical two-body (dotted line) and three-body (solid line) 
results, see also Ref. [5]. Concerning the shape of the en- 
ergy spectrum we have an excellent agreement between 
the experimental data and our results obtained using the 
three-body model, while the pure two-body calculation 
gives a too narrow peak. As to the absolute values, the 
three-body energy spectrum is approximately 20% above 
the experimental data. However, the most important 
lesson from this comparison is that for two different as- 
sumptions concerning the nuclear structure, but keeping 
the rms intercluster distance fixed, we obtain very differ- 
ent shapes of the calculated energy spectra. Thus, one 
can conclude that the interpretation of energy spectra is 
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highly model dependent. A final remark in connection to 
this observation is that, in order to interpret experimen- 
tal data correctly, it is very important to fix the spec- 
troscopic factors of different two-body and many-body 
components. Therefore, we want to stress the usefulness 
of experiments where EMD is studied in complete kine- 
matics. Examples of interesting channels in the 8 B case 
is 8 B -> 7 Be(l/2-) + p + 7 and 8 B -» 3 Hc + a + p. 
Some progress has already been made in this direction. 
Recently, Be fragments and 7-rays were measured in 
coincidence after breakup on a light target by Cortina- 
Gil et al. [11] which resulted in a clear observation of the 
excited core component of the WF. 




1.0 
E (MeV) 



FIG. 4: The 8 B EMD energy spectrum obtained at 
82.7 MeV/nucleon on Pb with 8 B scattering angles 9s B < 
1.77°. The data points are from [3]. The curves show our 
analytical two-body (dotted line) and three-body (solid line) 
results, and the numerical calculations: without FSI (long- 
dashed), and with FSI (short-dashed). All theoretical curves 
have been corrected for experimental resolution and accep- 
tance. 

We have also performed numerical calculations, based 
on first-order perturbation theory, where the numerical 
8 B bound state WF from [28] was used. In the first in- 
vestigation we assumed a pure Coulomb continuum and 
the results of this calculation is shown as a long-dashed 
line in Fig. 4. We note that the obtained energy spec- 
trum compares rather well with our analytical three-body 
model except for a 15% difference in the peak height. 
Since our analytical model has the correct asymptotic 
behavior for large r and a three-body behavior at small 
r, the main difference to the numerical WF should be 
in the intermediate region, and this is exactly the region 
which dominates the transition matrix elements for en- 
ergies corresponding to the peak of the energy spectrum. 
This fact explains the observed discrepancy. 

In our second numerical investigation we studied the 
influence of 7 Be + p FSI. As previously mentioned, the 
low-energy continuum of 8 B is dominated by positive par- 
ity states [25] which are only relevant for E2 (~10% of 
the total cross section) and Ml transitions. Note that the 
latter only plays a role in the vicinity of the narrow 1 + 
resonance at 0.64 MeV above threshold, and is therefore 



not included in our calculations. However, the possible 
existence of a very broad negative parity state at high 
excitation energy can still have a strong influence on the 
energy spectrum. Effects of such a state was observed 
in a recent elastic proton scattering experiment [27] from 
which the authors made a 2~ spin-parity assignment and, 
from an i?-matrix analysis, they obtained a best fit with 
the parameters E = 3.5 ± 0.5 MeV and T = 8 ± 4 MeV. 
We have included such a broad continuum structure by 
adding an attractive potential in the s-wave channel. The 
effects of this are clearly seen in Fig. 4 (short-dashed 
line): the total cross section is reduced and the contin- 
uum strength is redistributed towards smaller energies. 
We conclude this comparison by stating that a broad neg- 
ative parity structure in the high energy continuum has 
a non-negligible influence on the EMD energy spectrum 
for energies > 0.4 MeV and that the parameters of such 
a state are still to be determined with greater accuracy. 



B. Extraction of El contribution using the low-_E 
energy spectrum 

As we have seen in the previous section the main model 
uncertainties in the EMD analysis are: the asymptotic 
normalization constant which depends on (i) the radius 
in combination with (ii) the spectroscopic factors of dif- 
ferent two-body components; (iii) the underlying many- 
body structure, and finally; (iv) the FSI. Despite the 
difficulties connected with measuring the cross section at 
small relative energies, we still suggest to focus on the 
low--E part of the energy spectrum. In this way one can 
avoid uncertainties associated with FSI (unless there are 
resonances very close to threshold) and with the many- 
body behavior of the WF at small intercluster distances. 
The asymptotic normalization constant in combination 
with the spectroscopic factor will enter as an absolute 
normalization of the cross section. However, since this 
normalization affects all multipolc transitions equally it 
is possible to calculate the ratio of two different multi- 
poles with a very good precision. 

We have investigated the ratio 



R 



EMD 



(E) 



d<7EMp(El)/dg 
d<7 emd (El + E2) / dE ' 



(23) 



and found that it is almost model independent at small 
relative energies. To demonstrate this we will continue to 
use the EMD of 8 B (on a Pb target at 82.7 MeV/nucleon 
with 6 m ; n = 30 fm) as an example. Reusing the different 
model calculations from Sec. IV A we can investigate the 
sensitivity of the ratio to different model assumptions. 
Firstly, in Fig. 5 we can see that the difference between 
a two-body and a three-body approach is less than 3% 
in the region of interest. This low sensitivity can be ex- 
plained by the fact that the low-i? part of the spectrum 
mainly probes the large r asymptotics of the radial WF. 
Shown in Fig. 5 are also results from the numerical calcu- 
lations introduced in Sec. IV A. It is clearly seen that the 
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influence from FSI is almost negligible at small relative 
energies where the numerical results and our analytical 
three-body model seem to converge. Unfortunately, the 
numerical accuracy of our calculations becomes question- 
able at small energies and the calculated ratio is there- 
fore only plotted down to 0.1 McV. In this context we 
want to stress that in our analytical model we are able 
to calculate the relevant transition matrix elements for 
all energies, including the limit E — > 0. In contrast, the 
numerical approaches will run into problems for small en- 
ergies since the continuum WF will be extremely small 
at relevant interclustcr distances. 
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FIG. 5: Fraction of the calculated energy spectrum (EMD of 
82.7 MeV/nucleon 8 B on Pb with scattering angles 9&~Q < 
1.77°) attributed to El transitions, see Eq. (23). The curves 
show two-body (dotted line) and three-body (solid line) ana- 
lytical results while the dashed line shows results of numerical 
calculations including FSI. The inset shows the same curves 
for a broader energy range. 



Furthermore, as can be seen in Fig. 6, the sensitivity 
of the Remu(E) ratio (23) to the intercluster distance 
is also very small. This feature is also expected since 
the intercluster distance determines the asymptotic nor- 
malization of the WF which, in turn, cancels when the 
ratio is calculated. However, as can be seen in the insets 
of Figs. 5 and 6, the results obtained using the different 
models diverge with increasing energy. 

One final question is well-founded: Since the transition 
matrix elements for very small energies depend mainly on 
the tail of the bound state WF, is it still justified to use 
our model WF (13) which has merely an approximate 
description of the Whittaker tail? The result of our nu- 
merical calculation (remember that the numerical bound 
state WF has the correct asymptotics) presented in Fig. 5 
indicates that it is justified, since it seems to converge 
with the analytical model for small E. Furthermore, in 
the limit E — > we are actually able to solve the radial 
integral (10) exactly even for the "exact" WF (11) which 
has a Whittaker tail. We find that Remd(E — > 0) calcu- 
lated with the model WF and with the "exact" WF agree 
within 0.5%. This result gives an additional justification 
to the use of our model WF for calculating transition 
matrix elements at small E. 

In summary we have found that the calculated 
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FIG. 6: Fraction of the calculated energy spectrum (EMD of 
82.7 MeV/nucleon 8 B on Pb with scattering angles 9s B < 
1.77°) attributed to El transitions, see Eq. (23). The curves 
show two-body results characterized by two different inter- 
cluster distances: r rm s = 3.3 fm (dash-dotted) and r rms = 
5.4 fm (dotted). The inset shows the same curves for a broader 
energy range. 



-Remd(-E) ratio of the energy spectrum is almost model 
independent at small relative energies. In first-order per- 
turbation theory this ratio can be expressed as 



R 



EMD 



(E) 



dg E M D (El)/dg 
d ( 7 E MD(El + E2)/cLE 



n E 2(E. 



= 1 



1 + 



n E1 (E, 
n E i(E 7 )' n 




(24) 



where we have introduced the E2/E1 ratio of photo- 
dissociation cross sections 



T E2 



(E) 



rEl 



(25) 



Naturally, the ratio Remu(E) will depend on the exper- 
imental conditions such as beam energy and minimum 
impact parameter. However, this dependence enters only 
in the spectra of virtual photons which are easily calcu- 
lated for any experimental conditions {6 m i n , -Ebcam} (see, 
e.g., Eq. (4.7) of [23]). In contrast, the ratio of 8 B photo- 
dissociation cross sections r 1 {E) (25) does not depend on 
the experimental conditions. Therefore, we provide this 
ratio, calculated within our analytical three-body model, 
in Fig. 7. The curve shown in Fig. 7 can approximately 
be described by the formula 



_E2 



AE) 



(E) 



a. 



El 



(E) 



(-1.7290 + 3.4663e 15526iJ ) • 10" 



(26) 

which describes the calculated curve with an accuracy 
of < 0.5% in the region < E < 0.3 McV. From this 
formula it is easy to obtain i?EMD (E) using Eq. (24) 
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FIG. 7: Ratio of photo-dissociation cross sections calculated 
within our analytical three-body model, see Eq. (25). 



C. Extraction of El contribution using two 
different experimental conditions 



scattering angles (minimum impact parameter). The El 
contribution to one of the measurements can then be ob- 
tained with the formula 



dE 



' A W A < 2 ) " 

A7-( 2 )/p \™EMD wW/p \ fl(T EMD 



d# 



d£ 



A2) 



(2), 



r (i). 



(30) 



The advantage of this method is that information can 
be obtained directly from experimental data. However, 
it should be emphasized that the formula is only valid 
under the assumptions of first-order perturbation the- 
ory and straight-line trajectories. Thus, this method 
can only be used for EMD at relatively large beam en- 
ergies where events characterized by large impact pa- 
rameters have to be selected. Furthermore, the exper- 
imental conditions must be chosen so that the difference 
da^^/dE — d(7g 2 j^ D /di? is observable and larger than 
the experimental uncertainty. 



One of the beauties with the method of virtual pho- 
tons is the separation of reaction kinematics and nuclear 
excitation dynamics into the spectrum of virtual pho- 
tons and the photo-absorption cross section, respectively. 
This separation can be used as an alternative method to 
extract the El contribution from the measured cross sec- 
tion. The objective is to use the fact that the cross sec- 
tion will depend on beam energy and minimum impact 
parameter only through the spectra of virtual photons. 
First, let us introduce the notation 



nE\(E 7 , bmin, -Ebcam) 

E-y 



(27) 



where we have indicated that the spectra of virtual pho- 
tons are functions of the beam energy and minimum im- 
pact parameter. In reality this dependence enters in the 
adiabaticity parameter 



(28) 



Assuming that the total cross section is dominated by El 
and E2 transitions we find that the energy spectrum for 
given experimental conditions {b^ n , E^ am } is given by 



da 



(0 

EMD 



N£(E 7 )a?(E)+Ng>(E 7 )a™(E), (29) 



(i), 



dE 



where we would like to remind the reader of the relation 
E = E 1 — Eq. Now we can use the fact that the vir- 
tual photon spectra depends differently on £ for different 
multipoles to extract the contribution from one of the 
multipoles. Let us assume that we have two sets of ex- 
perimental data from the same experimental setup; the 
only difference being the beam energy and/or the selected 



V. CONCLUSION 

In this paper we have studied electromagnetic pro- 
cesses involving loosely bound nuclei. To this aim we 
have developed an analytical model which is based on 
the use of radial model functions that give a realistic de- 
scription of two-body WFs (or the two-body projections 
of many-body WFs) for all radii, see also [5]. We have 
used this model to study EMD of 8 B, but have also indi- 
cated how it can be applied to other reactions and other 
nuclei. For example, it should provide an important tool 
to investigate the low-energy behavior of the astrophysi- 
cal S- factor. 

We have also presented numerical calculations based 
on the three-body model of 8 B, developed in Ref. [28], 
and on recent experimental information on a broad neg- 
ative parity state in the 8 B continuum. Combining the 
results of our analytical model, and of these numerical 
calculations, has allowed us to investigate the sensitivity 
of calculated EMD energy spectra to different model as- 
sumptions. We concluded that the magnitude of the cross 
section depends strongly on the intercluster distance and 
on the spectroscopic factors of different many-body chan- 
nels; while the shape of the energy spectrum is very sen- 
sitive to the few-body structure. Finally, we found that a 
broad negative parity state at high excitation energy will 
influence the energy spectrum for energies > 0.4 MeV, 
and will lead to a reduction of the cross section. 

However, the main purpose of this paper has been to 
investigate the problem of how to extract the El con- 
tribution from a measured EMD energy spectrum. This 
question is of great significance for the gathering of infor- 
mation on astrophysically interesting radiative capture 
reactions from EMD experiments (note the relation be- 
tween the El component of the EMD energy spectrum 
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and the astrophysical S- factor via Eqs. (4-6)). The main 
method so far has been to study the asymmetries in angu- 
lar or momentum distributions. However, this asymme- 
try (which is due to E1-E2 interference) depends strongly 
on details of the FSI which, in turn, are often relatively 
unknown. Furthermore, the E1-E2 interference terms do 
not themselves contribute to the integrated cross sections 
to which the S'-factor is related. Instead, we have pro- 
posed two novel, and less model dependent, approaches 
to extract the El contribution from a measured EMD 
energy spectrum: (i) Firstly, we demonstrated that the 
ratio of EMD cross sections ctemd(E1)/ctemd(E1 + E2) is 
almost model independent at small relative energies. We 
also provided an analytical formula to calculate this ra- 
tio for any experimental conditions, (ii) Secondly, we 
demonstrated how two sets of experimental data, ob- 
tained with different -Eboam and/or 6 m i n , can be used to 
extract the El component. This method relies on the 
fact that the strengths of different multipole components 



depend on the beam energy and minimum impact pa- 
rameter, and in first-order perturbation theory this de- 
pendence enters only in the virtual photon spectra. 

Since the proposed two methods are not directly con- 
nected to each other, they can be used independently and 
the results can be compared to each other. However, both 
methods, but in particular the first one, require that the 
energy spectrum is measured down to very small relative 
energies (100-300 keV), which will probably prove to be 
a difficult challenge. 
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